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ABSTRACT 


Hll regions are particularly interesting because they can generate dense layers 
of gas and dust, elongated columns or pillars of gas pointing towards the ioniz¬ 
ing sources, and cometary globules of dense gas, where triggered star formation 
can occur. Understanding the interplay between the ionizing radiation and the 
dense surrounding gas is very important to explain the origin of these peculiar 
structures, and hence to characterize triggered star formation. G46.5-0.2 (G46), 
a poorly studied galactic HII region located at about 4 kpc, is an excellent target 
to perform this kind of studies. Using public molecular data extracted from the 
Galactic Ring Survey ( 13 CO J=l-0) and from the James Clerk Maxwell Telescope 
data archive ( 12 CO, 13 CO, C 18 0 J=3-2, HCO + and HCN J=4-3), and infrared 
data from the GLIMPSE and MIPSGAL surveys, we perform a complete study 
of G46, its molecular environment and the young stellar objects placed around 
it. We found that G46, probably excited by an 07V star, is located close to 
the edge of the GRSMC G046.34-00.21 molecular cloud. It presents a horse-shoe 
morphology opening in direction of the cloud. We observed a filamentary struc¬ 
ture in the molecular gas likely related to G46 and not considerable molecular 
emission towards its open border. We found that about 10' towards the southwest 
of G46 there are some pillar-like features, shining at 8 //in and pointing towards 
the Hll region open border. We propose that the pillar-like features were carved 
and sculpted by the ionizing flux from G46. We found several young stellar ob¬ 
jects likely embedded in the molecular cloud grouped in two main concentrations: 
one, closer to the G46 open border consisting of Class II type sources, and other 
one mostly composed by Class I type YSOs located just ahead the pillars-like 
features, strongly suggesting an age gradient in the YSOs distribution. 



Subject headings: ISM: Hu regions - individual: G46.5-0.2 - ISM: molecular clouds 
individual: GRSMC G046.34-00.21 - Stars: formation ISM: molecules 
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Introduction 


Massive stars play a key role in the evolution of the Galaxy. They are the principal 
source of heavy elements and UV radiation. Through the combination of winds, massive 
outflows, expanding Hll regions and supernova explosions, they provide an important 
source of mixing turbulence in the interstellar medium (ISM). Hll regions are particularly 
interesting because they can trig ger t he formation of a new generation of stars. Processes 


such ‘collect and co llapse’ (e.g. 


(e.g. 


Bertoldi 


Elme gre en et ah 


1995) and radiative driven implosion (RDI) 


19891) have been con vincingly demonstrated that occur around Hu region s 


leading to formation of st 

ars (e.g. 

Ortega et al. 

2013; 

Dirienzo et al. 

2012- 

2011; 

Zavagno et al. 

2010; 

Deharveng et al. 

2008 

ThoniDSon et al. 

2004a 

)• 


Paron et al. 


generate dense layers of gas and dust at the interface between the Hu regions and their 
parent molecular cloud, elongated columns (pill ars) o f gas p ointing towards the ionizing 


sources, and cometary globules of dense gas (Trcmblin et al. 


20131) . Understanding the 


interplay between the ionizing radiation and the dense surrounding gas is very important 
to explain the origin of these peculiar structures, and hence to characterize triggered star 
formation. 

G46.5-0.2 (hereafter G46) is a poorly studied Hn region about 8' in size located at RA 


= 19 h 17 m 26fi, dec = 11° 55(5T( (J2000), on the b order of t 


lie mo 


GRSMC G046.34-00.21 ( Rathborne et al. 


20091 ) 


Lockman 


ecular cloud catalogued as 


(19891) ob serve d a recombination 


line at v LSR rs j 57.2 km s 1 towards this region, 
absorption features up to 56.6 km s' 1 and based on 
for this Hll region the near distance of 3.8 kpc, while 


Kuchar & Bania ( 19941 ) detected HI 


dnematical considerations proposed 
, using the same 


Quireza et al 


techniques, favored the far distance of 7.8 kpc. Finally, 


(2006) 


Anderson & Bania ( 2009 ) resolved 


the kinematic distance ambiguity for G46 based on existing HI and 13 CO sky surveys, 
confirming a near distance of 4 kpc, which will be adopted in what follows. 
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Figure [I] shows a two-color composite image of G46 field where the Spitzer -IRAC 8 
pm emission distribution is displayed in cyan and the Spzfzer-MIPSGAL at 24 pm in red. 
The emission displayed in cyan mainly corresponds to radiation originated in polycyclic 
aromatic hydrocarbon (PAH) molecules. The far-UV photons leaking from the Hll region 
excite the PAHs, which then emit in the mid-infrared. This emission encircles the 24 pm 
emission (in red) that arises from very small size dust grains which are heated but not 
destroyed by the UV photons. The emission at 8 pm has a horse-shoe morphology, being 
thicker and more intense towards the east-northeast direction and open to the southwest, 
where it dilutes into the surrounding medium. Precisely in the southwest direction, about 
10' away from G46 there are pillar-like features shining at 8 pm. In this paper we explore 
the interrelationship among the Hll region G46, the environmental conditions and the 
young stellar objects in the region, underscoring evidence of triggered stellar formation in 
the immediate vicinity of G46 and farther. 


Data 


2.1. Molecular data 


The 13 CO J=l—0 data were extracted from the Galactic Ring Survey (GRS). The 
survey, performed by the Boston University and the Five College Radio Astronomy 
Observatory, maps the Galactic Ring in the mentioned molecular line with an angular 


and spectral resolutions of 46" and 0.2 km s 1 , respectively (see 


Jackson et al. 


2006 b The 


observations were performed in both position-switching and on-the-fly mapping modes, 
achieving an angular sampling of 22". Data are presented in main beam temperature ( T m b ) 


Additionally we used 12 CO, 13 CO, C 18 0 J=3 2, and HCO + , HCN J=4-3 data extracted 





from the James Clerk Maxwell Telescope (JCMT) data archive^]. These observations 
(Proposal I.D M10AH02, PJ. Jonathan Williams) were carried out with the JCMT in 
Mauna Kea, Hawaii, using the HARP-ACSIS instrument. The angular and spectral 
resolutions are about 15" and 0.05 km s -1 . We used the reduced data (cal. level 2), which 
were reduced using the standard ORAC-DR pipelines 1 2 3 ]. Data are presented in units of 
corrected antenna temperature T*, which is related to the main beam temperature (T mb ) 


using T mb = TX/rj mb . Following 


Buckle et al. 


(2009) a value of rj mb = 0.6 was used. The 


zenith opacity was between 0.04 and 0.06 for all the observed lines. Even though we 
extracted reduced data, in the case of C 18 0 J=3 2, HCO + , and HCN J=4-3 data we 
applied a Hanning smoothing to improve the signal-to-noise ratio which altered the spectral 
resolution in no more than a factor of 2, and minor polynomials were used for baseline 


corrections using XSpec 




2.2. Infrared data 


In this work, archived infrared data are also used to reveal star formation in G46. 
These data come from the GLIMPSE and MIPSGAL surveys. 


GLIMPSE (Galactic Legacy Infared Mid-Plane Survey Extraodbiaire. 


20031) using the IRAC ( using the InfraRed 


Spitzer Space Telescope flWerner et al 


Array Camera, 


Fazio et al 


Jen i am in et al 


2004) on board the 


2004) surveyed the inner 130 degrees of the Galactic 


Plane at 3.6, 4.5, 5.8, and 8.0 pm bands. The 5 cr sensitivities of the four bands are 0.2, 0.2, 


1 www. j ach. hawaii. edu/J CM T /archive / 

2 http: / / www. j ach. hawaii .edu /JCMT/archive / C AD C_quickguide. lit ml 

3 XSpec is a spectral line reduction package for astronomy which has been developed by 


Per Bergman at Onsala Space Observatory. 
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0.4, 0.4 mJy, respectively. In addition to images, GLIMPSE survey performed point-source 


photometry. Photometric data at the J (1.25 /ini), H (1.65 /mi), an d Ks ( 2.17 pm) 


lands 


from the Two Micron All Sky Point Source Catalog (2MASS PSC, 


Skrutskie et ah 


2006) 


are provided in the GLIMPSE Catalog to build-up a 7-bands photometric system. Both the 
images and the catalog of GLIMPSE are publicly available at the InfraRed Science Archive 
(irsa£ where we have retrieved cutouts in the four IRAC bands and point sources in 
a 31' x 24' region centered at 0:2000 = 19 h 17 m 18 s .7, h 2 ooo = +11°52'22".9. We restricted 
the extracted catalog to be a more reliable data set with the following criteria: a) only 
sources with photometric errors no larger than 0.2 mag in the IRAC bands are taken into 
account; b) for 2MASS bands, a threshold of 0.1 mag photometric error is used to forsake 
unreliable photometric val ues. The photometric errors for IRAC and 2MASS sources have 


been adopted following 


Gutermuth et al. 


(120091 1. 


MIPSGAL (Carey et al. 


2009) is a complement to the GLIMPSE le gacy su rvey. 


survey using the MIPS (Multiband Infrared Photometer for Spitzer, 


Rieke et al. 


2004) 


'his 


instrument on board the Spitzer Space Telescope surveyed an area comparable to that of 
GLIMPSE. The version 3.0 of MIPSGAL data includes mosaics only at 24 pm with sky 
coverage of |6| < 1° for —68°< l < 69°, and |6| < 3° for —8°< l < 9°. The spatial resolution 
and 5cr sensitivity at 24 pm are 6" and 1.7 mJy. From the IRS A server, we have extracted 
a cutout of a region same as that of GLIMPSE cutouts. 


We conducted point-source extraction and aperture photometry of point sources, in 


the 


MIPS 24 pm image using the PSF fitting capability of IRAF/DAOPHOT ( Stetso n 


1987 ). The PSF was determined to be about 6" by fitting the profiles of ten bright point 


sources in the investigated field. The standard deviation (er) of the sky was estimated to 
be about 5.0 x 10 -6 Jy pixel -1 , daofind was used to extract candidates of point sources 


4 http://irsa. ipac.caltech.edu/index.html 



























with a threshold of 10 a. The final sources were determined by visual inspection using the 
tvmark task. Sources affected by ghosts, diffraction spikes, halos from bright sources and 
artifacts residing in bright extended emission were rejected. For aperture photometry of the 
extracted sources, radii of the apertures and inner and outer limits of the sky annuli were 
selected to be 4.5", 15" and 21", respectively. Magnitudes of the extracted sources were 
determined using the magnitude zero point of 7.17 Jy provided in the MIPS instrument 
handbook. Finally, we cross-matched the 24 /mi sources with the GLIMPSE catalog using 
a cone radius of 2". 


Additionally, we have queried the MIPSGAL point source catalog given by 


Gutermuth & Heyer 020151) . A total of 309 sources are found in the region of interest. 


Among them, 305 have been covered by our 24 /mi point source list. The other four sources 
are weaker than 11.5 mag. A comparison of the fluxes indicates that our photometry is 


consistent with that given in 
than 5% in magnitude. 


Gutermuth & Heyer 


(2015), the mean difference is smaller 


The completeness of our catalog was estimated by counting the number of sources as 
a function of magnitude. A histogram plot of source magnitudes at each band has been 
plotted and carefully inspected. The sources magnitudes are exponentially distributed. We 
plot the number of sources in log scale and the magnitude in linear scale, then a straight 
line can be fitted. The magnitude at which a deviation emerges has been considered as the 
completeness limit. We found that our catalog is complete to a magnitude of 14.0 at 3.6 
H m, 13.5 at 4.5 /tm, 12.3 at 5.8 /mi, 12.0 at 8.0 /mi, and 8.5 at 24 /mi. For the central 
region with diffuse mid-infrared emission, our catalog is complete to a magnitude of 14.0 at 
3.6 /mi, 13.5 at 4.5 /mi, 12.0 at 5.8 /mi, 11.5 at 8.0 /mi, and 8.0 at 24 /mi. 


5 http://irsa. ipac.caltech.edu/data/SPITZER/docs/niips/ 
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3. Results and discussion 

In what follows we separately analyze the physical properties of the Hll region G46 
(Sec. 3.1), of the ambient molecular gas (Sect. 3.2) and of the young stellar objects (YSOs) 
near G46 (Sec. 3.3), with the aim of investigating the impact of the ionizing radiation on 
the surrounding medium and the likelihood of triggered star formation. 


3.1. Exciting star(s) candidate(s) and morphology of G46 


able OB-type stars catalogues 


A search for the 

exciting star of G46 in the 

avai 

Mafz Apcllaniz et al. 

2013; 

Mafz-Apellaniz et al. 

2004 


Reed 


2003)), was unsuccessful. We 


can, anyway, conjecture the earliest spectral type of a probable single ionizing star from the 
radio continuum emission of the Hll region. The number of UV ionizing ph otons nee ded t o 
keep an Hll region ionized is given by N uv = 0.76 x 10 47 T 4 ~°' 45 S v D/ pc ( Chaisson 1976 1. 
where T 4 is the electron temperature in units of 10 4 K, D kpc the distance in kpc, z/ G Hz 
the frequency in GHz, and S„ the measured total flux density in Jy. Assuming a typical 


electron temperature of T = 1C 


4 K, a distance of 4 kpc, and a total flux density of 2.7 Jy 


as measured from the MAGPISfJ image at 1.4 GHz, the total amount of ionizing photons 
needed to keep this source ionized turns out to be about N w = (3.5 ± 1.2) x 10 48 ph s _1 . 
Assuming errors of about 10% in both the distance and the radio continuum flux density, 


we conclude that the exciting star should be an 07V ({Martins et al 


2005) 


In a rough attempt to identify the star(s) candidate(s) responsible for the ionized gas in 


the region we performed an optical and infrared photometric study of the qoin 


the area (see Fig. [2]) based on the astrometric UCAC3 Catalog (IZacharias et ah 


sources m 


20101 1. Only 


sources with detection in the optical B and R bands, and in the three near-infrared (NIR) 


6 http://third.ucllnl.org/gps/ 
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J,H,K bands extracted from the Two Micron All Sky Survey (2MASS)Q, were considered. 
We found eight sources inside a circle of size 4' centered at the central coordinate of the 
Hll region. Their locations are shown in Fig. [21 From the optical and infrared magnitudes 
we constructed their spectral energy distribution (SED). We fit the available magnitudes 


(optical B and R bands, the three JHK 2MASS banc 


using the Kurucz phot ospheric models (IKurucz 


s, and the four Spitzer-IRAC bands) 


Robitaille et ah 


197911 included in the tool developed by 


(l2007ir l to obtain the effective temperature, T e ff, of each source. The 


fitting tool requires the assumption of the visual extinction, A„, and the distance. We 
adopt a distance between 3.5 and 4.5 kpc. Regarding the A v , we derived it for each 
source from their (J—H) and (H—K) colors. We assumed the interstellar reddening law of 
Rieke fe Lebofskv ( 1985 1 (A j /Ay = 0.282, A h /Ay =0.175, and A k /Ay = 0.112) and 
the intrinsic colors (J—H ) 0 and (H—K ) 0 obtained by Mart i ns fe Plez (12006 b 


Among the eight stars found in the region, only three are compatible with a massive 
star located at the distance of G46. The effective temperatures of the massive candidates 
stars obtained from the fitting are shown in Table [H Stars #6 and ^8 have a T e ff of 
about 30000 K and 35000 K, respect ively, which agree with the temperatures of 09.5V 


and 07V stars, respectively (ISchaerer fe de Koter 


1997( 1. The models predict an effective 


temperature of ~ 26000 K for source #3, suggesting that this star would be of a spectral 
type later than B0. The fitted SEDs for sources ^6 and ^8 are shown in Fig.[3] For the 
remaining sources #1, #2, #4, #5 and #7, the Kurucz’s models fail to fit with confidence 
their optical and NIR magnitudes as a massive star at about 4 kpc. These sources are 


'2MASS is a joint project of the University of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Technology, funded by the National Aeronautics 
and Space Administration and the National Science Foundation. 

8 http://caravan.astro, wisc.edu/protostars/ 
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probably less massive foreground stars not related with G46. Thus, stars 7^6 and 7^8 are the 
most probable candidates to be the exciting stars of the Hu region. Finally, the location of 
the star 7^8 with respect to the radio continuum emission (red in Fig. [2]), which looks like 
an incomplete ring, seems to better explain the curved morphology observed in the ionized 
gas distribution. 


3.2. Molecular gas analysis 


The Hu region G46 is located in 
GRSMC G046.34-00.21 catal ogued by 


54.3 ± 3.9 km s 


-1 


projection onto a bord er of the molecular cloud 


Rathborne et al. 


Roman-Duval et a! 


(120091 ) with an associated v LSR of 


( 20091 ) established a kinematical distance of about 


4 kpc for this molecular cloud, the same as for the fill region G46, supporting that both 
objects are linked. Figure|4] (left) displays the 13 CO J=l-0 integrated emission distribution 
in the 40 - 62 km s -1 velocity range. It can be seen that the Hll region seems to be evolving 
on the northeastern edge of the cloud. The molecular gas appears located in projection 
onto the diffuse emission at 8 /irn observed towards the southwest of the open border of 
G46. It is noteworthy that there is no evidence of related molecular emission towards the 
east-northeast edge of the Hii region, where its associated PAHs emission is thicker and 
more intense. 


The numerous studies of Hll regions and their surroundings usually show that the 
associated PAHs emission is more intense towards the borders of the Hll regions with 
presence of molecular material. I 11 other words, it is expected that an ionized bubble opens 
in the direction away from the cloud edge because the ionized gas would be less confined 
and could st ream-out in t o the lower density ISM, forming what is known as a “blister-type” 


Hii region ([Isra el 


1978; 


Tenorio-Tagle 


1979). In the G46 case, however, the ring opens to 


the southwest, facing the elongated molecular cloud where apparently the density is higher. 
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A possible explanation can be found from the morphology of the G46 radio continuum 
emission, which suggests a stalling effect of the ionized gas against the photodissociation 
region (PDR) towards the northeast (see Fig. [2]). This fact suggests that although there 
is no 13 CO J=l—0 detected towards this region, the density of the photodissociated gas is 
enough to confine the ionized gas. The atypical thickness of the PDR towards this edge 
could be due to the presence of relatively low density gas which is expected towards the 
edge of a molecular cloud whereby the far-UV photons would have greater penetration in 
the material. By the other hand, the non detection of radio continuum emission towards 
the open border of G46 suggests that the ionized gas escapes from G46 and dilutes in the 
ISM. The presence of diffuse emission at 8 /jrn towards the southwest of G46 supports this 
scenario. 


3.2.1. The pillar-like features 


As above mentioned, about 10' towards the southwest of the Hii region (~12 pc at the 
assumed distance of 4 kpc) it can be seen some pillar-like features emitting at 8 /jrn, which 
are pointing to the ionizing source of G46 and seem to be embedded in the molecular cloud 
GRSMC G046.34-00.21. Even though it seems to be a large distance, several pillars were 
found_as far from the ionizing source as those analyzed in this work (e.g in the Vulpecula 


rift; 


Billot et ah 


2010). What it is really remarkable in this case is that these pillars-like 


features, resembling heaps and corrugations, are well outside of the Hii region, in contrary 
as is usually found: bubble-like structures with pillars inside or over their boundaries. 
Taking into account that the pillars point to the G46 open border, it is suggested that they 
were produced by UV photons escaped from the Hu region. 


By inspecting the 13 CO J=l-0 towards the pillar-like features (Fig. |4}- right), we find 
that there are two molecular structures associated with them. The easternmost IR pillar-like 
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feature has associated molecular gas between 54.0 and 56.5 km s^ 1 (white contours in 
Fig.HJ-right), while the western one is related to a molecular structure ranging between 56.5 
and 59.0 km s -1 (yellow contours in Fig. fright). The different velocity ranges for these 
molecular structures may either indicate that both features are located at slightly different 
positions along the line of sight or that they have a somewhat different kinematics. It 
is remarkable the good agreement between the morphology of the molecular gas and the 
pillar-like features as seen in IR. Moreover, the molecular structures present a morphology 
consisting of a dense head with a less dense tail as is usually found and predicted by 


models and o 


oserva 


(Pound et al. 


2007 


ions towards this k i nd of structures in the surroundings of Hu regions 


Schuller et al. 


2006 


Pound et al. 


2005)- However, following the work of 


Macke y & Lim (1201011 . the appearance of these structures resembles a previous evolutive 


stage in the formation of pillars because it can not be appreciated a well-formed tail behind 
their heads. In what follow we characterize the interaction between the ionization radiation 
escaping from G46 and the pillar-like features. 

In order to study the radiation influence over the tips of the pillar-like features and the 


possibilities of triggered 


star formation via_radiative_driven implosion (RDI, e.g. 


1989, 

Lefloch et al. 

1997, 

Kessel-Devnet & Burkert 

CO 

o 

o 

CN 


Bertoldi 


2003 ) we evaluate the pressure balance 


between the ionized gas stalling at the head of the pillars and the neutral gas of their 
interiors. Assuming that the ionizing photons came from an 07V-type star located at 
the center of G46 (see S ect .137 


star from 


Schaerer fc de Koter 


, we use the predicted ionizing photon flux for an 07V 


( 19971 ) and the projected distance between the star and 


the pillars to roughly estimate the amount of UV photons arriving at the surface of the 
pillars in 4> pre ~ 2 x 10 8 cm -2 s -1 . This value represents an upper limit due to the fact 
that the projected distance between the star and the pillars is a lower limit to the actual 
distance between them. Our predicted ionizing photon flux is similar to those measured 


towards several bright-rimmed clouds ({Thompson et al. 


2004a.b), which like the pillars, are 
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molecular clouds sculpted by the radiation leaking from HlI regions, and in many cases they 


have the same large scale morphology as the pillars (Thompson et al 


2004a) 


Using the above obtained & pre and following 


Thompson et al.l (2 004b ) we estimate an 


upper limit for the electron density of about 64 cm 3 for the expected ionized boundary 
layer (IBL) at the tip of both pillar-like structures. This is almost 3 times greater than 


the critical value o 


(Lefloch & Lazaret! 


25 cm 3 above which an IBL is able to develop around a cloud 


1994) 1. To obtain this value we consider that the pillar heads have a 


radius of 0.5' and assume an effective thickness of the ionized boundary layer of t] = 0.1. 
Then, using a typical sound speed of the ionized gas of 11.4 km s _1 , we obtain that the 
pillars tips are supporting an external pressure of P ex t/k ~ 2.1 x 10 6 cm -3 K. It is important 


to note that the obtained ~P ext represents strictly an upper 


an upper limit and due to the used r) (see 


Thompson el al. 


imit because the predicted <f> is 


2004b) 


On the other side, integrating the 13 CO J=l-0 emission over the area that contain 
each pillar head and using the typical LTE formulae to derive the 13 CO and Lp column 


den sities by assuming a T ex = 10 K and [H 2 / 13 CO] = 5 x 10 5 (see e.g. 


Yamaguchi et al. 


19991 1. we estimate the mass of each pillar head in about 450 and 300 M 0 for the eastern 
and western one, respectively. The area of each pillar head was assumed by taking into 
account the radius of their tips curvature, this is a circular area with radius of 0.5' for both 
pillar heads. Thus, assuming spherical shapes, we estimate the molecular densities in about 
7.0 x 10 3 and 4.6 x 10 3 cm -3 , respectively. Then, considering the velocity interval in which 
each pillar structure extends, we obtain the velocity dispersion a v ~ 1.05_km s -1 , wh ich 


is in 

agreement with 

2010 ; 

Dale et al. 

2012 ) 


hose predicted by models of pillars formation (jGr itschneder et ah 


pressure for each pillar head in P i V ±/k ~ 1.8 x 1 0, and ~ 1.2 x 10 6 cm 3 K for the eastern 


and western one, respectively. As 


Thompson et al.l (j2004bf) state, these pressure values are 
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very likely underestimated because the 13 CO J=l-0 line underestimates the true density 
and this molecule may be depleted by selective photodissociation at the boundary of the 
clouds. Following these authors, the internal pressure is likely underestimated by no more 
than a factor of 15, thus we conclude that P i n t/k should be between 2 x 10 6 and 3 x 10' 
cm -3 K. 


In conclusion, we obtain P int > P ext , which suggests that the diluted ionization front 
stalls at the pillar heads, probably until the effects of mass evaporation and increasing 


recombination within the IB 


pressure (Lefloch & Lazarefl 


raises the ionized gas pressure to equilibrium with the interior 


1994). This result shows that it is unlikely that a shock is 


propagating farther into the molecular gas, discarding that the RDI mechanism is on going 
in the pillars interior. 


3.2.2. The molecular gas towards GJ^6 open border 

With the purpose of investigating the molecular material towards the open border of 
G46 we used higher-angular resolution data of several molecular species extracted from the 
JCMT database acquired towards this region (see the rectangle in Fig. 01-left). Figure [5] 
shows the integrated emission of the 12 CO, 13 CO, and C 18 0 J=3-2, while Fig. [6] displays 
the HCO + and HCN J=4-3 integrated emission. The 12 CO and 13 CO J=3-2 emission 
show clumpy, elongated molecular features, with a high density filament extending towards 
the southwest of the mapped region. Towards the open border of G46 do not appear any 
considerable molecular emission, suggesting the presence of a pre-existing region with scarce 
molecular gas or that the UV-photons have carved the molecular cloud. In any case this can 
be the path followed by the UV-photons escaped from G46 to reach the farther pillar-like 


structures. 
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The C 18 0 J=3-2, HCO+ and HCN J=4-3 emissions are concentrated in a small 


compact clump related to the infrared dark cloud IRDC 046.424-0.237 (Peretto & Fuller 


20091 ) and to the millimet er continuum source BGPS G046.427-00.237 (see the Bolocam 


Galactic Plane Survey v2; 


Gi nsbu rg et ah 


20 13 1). suggesting that the emission lines emanate 


from a deeply embedded clump. Moreover, this structure is associated with the cold 


high-mass clump G46.43-0.24 flWienen et al. 


2012 ) in which the (1,1) and (2,2) ammonia 


lines were detected at vlsr ~ 52.3 and 52.9 km s 1 , respectively. The detection of these 
molecular species, mainly the HCN J=4-3 line that its criti cal density can be between some 


10° and 10 s cm 3 (iTakakuwa et al 


2007|; 


Greve et al 


2009), indicates the presence of very 


high density gas. Taking into account that this clump is a potential site of star formation, 
in what follows we analyze it. 


The line parameters of the observed molecular transitions towards the center of the 
clump are given in Tablc[2]as derived from Gaussian fits from the spectra shown in Fig. [7] 
The 12 CO J=3-2 was fitted using three Gaussians, while the parameters of the others lines 
were obtained from single-components Gaussian fits which coincide in velocity with the 
main 12 CO component. The different velocity components observed in the 12 CO spectrum 
are present almost in the whole region. Thus, we conclude that they correspond to different 
molecular components seen along the line of sight, reflecting the dumpiness in the region. 


In order to have a rough estimate of the molecular clump mass we assume local 
thermodynamic equilibrium (LTE). We calculate the excitation temperature from 


Tex (3 —> 2) — 


16.59K 


( 1 ) 


where T„ . 

obtaining T £ 


ln[l + 16.59K/(T max ( 12 CO) + 0.036K)] 

( 12 CO) is the 12 CO peak temperature towards the clump center at ~ 52 km s -1 
18 K. We derive the 13 CO and C 18 0 optical depths r i3 and Ti 8 from (e.g. 


Buckle et al. 


2010 ): 


13 T mfe _ 1 - exp(-r 13 ) 

18 T mb 1 - exp(-T 13 /X) ’ 


( 2 ) 
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where 13 T m b and 18 T, n b are the peak temperatures of the 13 C0 and C 18 0 J=3-2 line at the 
center of the region, and X = 8.4 is the assumed isotope abundance ratio [ 13 C0/C 18 0] 


(Frerking et ah 


1982 


Wilson 


19991), obtaining r i3 ~ 3.5 and r 18 ~ 0.4, which indicate that 


the C 18 0 J=3-2 line appears moderately optically thin. Thus, we estimate its column 
density from: 


N(C 18 0) = 8.26 x 10 13 e 


13 


i5.8i T ex + 0.88 


1 — e Te; 


J{T ex ) - J (Tbg) 


- mb 


dv 


( 3 ) 


with 


J(T) = 


hv/k 


( 4 ) 


To obtain the molecular 
[H 2 /C 18 0] = 5.88 x 10 6 (Frerking et al. 
from: 


ex p(w ) - 1 

lydrogen column density N(H 2 ) w e assume an abundance ratio of 


1982 


Wilson 19991). Finally the mass was derived 


M — n m H [D 2 Qi Nj(H 2 )], (5) 

i 

where 11 is the solid angle subtended by the beam size, rri\\ is the hydrogen mass, /i, the 
mean molecular weight, is assumed to be 2.8 by taking into account a relative helium 
abundance of 25 %, and D is the distance. Summation was performed over all beam 
positions on the molecular structure observed in C 18 0 displayed in contours in Fig. [5] The 
obtained mass is about 375 M©. Assuming an ellipsoidal volume with semi-axis of 45", 
20", and 20" we derive a density of n ~ 1.5 x 10 4 cm -3 . The density should be higher in 
the innermost region of the clump where the HCN J=4-3 line emanates, showing a density 
gradient. Using the deconvolved radius of about 0.5 pc calculated from: 


R, 


clump 


S — beam area 


7r 


( 6 ) 


where S is the area inside the clump, we appreciate that the above obtained mass value 
is slightly ove r the mass-s ize threshold for massive star formation in IRDCs presented in 
Ka.uffma.nn & Pillai (2010), suggesting that massive YSOs can be formed within this clump. 
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Using the derived HCO + and HCN J=4-3 parameters listed in Tablc[2]we performed 
a non-LTE study of these molecular species with the code RADEX, whic h uses the mean 


escap e probability approximation for the radiative transfer equation (van der Tak et al 


2007). Using the measured Av we ran th e code to fid T m b and estimate the column 


Wienen et al. 


( 20121 ) measured from the ammonia lines 


densities. Taking into account that 
a kinetic temperature of ~ 16 K towards this clump, we fix this parameter in 20 K 
and assume densities between 10 5 and 10' cm -3 to obtain the values presented in Table[3j 
From the obtained column densities we observe abundance HCN/HCCU ratios of about 
4.4, 3.8, and 1.2 for iih 2 = 10 5 , 10 6 , and 10' cm -3 , respectively. In all cases we obtain a 


NfHCNf/NfHCO+f ratio 


G48.66-0.22 (Pitann et al. 


arger than unity as found towards several clumps in the IRDC 


20131) and the active star-forming region W49A ( R o berts et al, 


2011 b As these authors point out, the steady-state chemical models for molecular species 


in gas-phase predict HCN > HCO + only for < 25 K, with a density of iih 2 = 10 6 cm 3 , 
which is consistent with our results. 


3.3. Star formation around G46 

Finally, given that the scenario is very favorable to initiate the formation of new 
generations of stars, we analyze in this section the existence and properties of young stellar 
objects in all the investigated area. 


3.3.1. Identification of YSOs 

Young stellar objects (YSOs) always show an excess in the infrared emission. The level 
of excess in the infrared can be effectively used for discriminating YSOs from field stars 
and distinguishing different evolutionary stages. At a considerably early evolutionary stage, 
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protostars are mostly embedded in dust envelopes, they exhibit large excess of infrared 


emission and an infrared spectral index cpR > —0.3 indicative of fla 


energy distribution (SED) at wavelength longward of 2 /mi flLada 


or ascending spec 


1987|; 


Greene et a! 


ral 


1994 1. For pre-main sequence (PMS) stars which possess optically thick disks, the SEDs 
tend to descend and the infrared spectral indices are in the range —1.6 < cpR < —0.3. 
Finally “transition disk” sources are more evolved YSOs, where the inner part of the disks 
have been cleared by photoevaporation of the central stars or by planet forming processes. 
For such YSOs it is expected to detect 1R excess at wavelengths longer than 16 /mi 


fjStrom et ah 


1989 ). These properties of YSOs make photometric observations in the near- 


to mid-infrared plausible to discriminate them from held stars. Different color-based source 
identification and classification schemes have been developed and verified in the practice. In 


this section, we identify potential YSOs following the scheme proposed by 


Gutermuth et al. 


020 091 ). The resulting YSOs are classified into Class I (protostars, including Class 0, Class 
I, and “hat spectrum” sources), Class II, and “transition disk” sources. 

There are several kinds of contaminants that could be misidentihed as YSOs in our 
original sample. Extragalactic contaminations could stem from star-forming galaxies and 


broad-line active galactic nuclei (AGN] 


red 5.8 and 8.0 /ini colors ([Stem et ah 


, whi ch sh ow PA H feature em ission yielding very 


2005 


Gutermuth et ah 


20081) . In our own Galaxy, 


unresolved knots of shock emission and resolved PAH emission are often detected in the 


IRAC ban ds, yieldin g ad c 


m 


Gutermuth et ah 


itional contaminations. Following the techniques presented 


( 2009 ) we exclude contaminants, and then identify and categorize 


potential YSOs as described as follows. 


In the hrst phase, only sources with valid detections in all four IRAC bands have been 
considered. Any source fulfilling color criteria of [3.6] — [4.5] > 0.7 and [4.5] — [5.8] > 0.7 
is regarded as a Class I YSO. In the remaining pool, Class II sources have been picked out 
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based on the constraints of i) [3.6] — [4.5] — o\ > 0.15, ii) [3.6] — [5.8] — <r 2 > 0.35, in) 
[4.5] - [8.0] - a 3 > 0.5, and m) [3.6] - [5.8] + a 2 < §df x (([4.5] - [8.0] - a 3 ) - 0.5) + 0.5 
Here, ay = cr([3.6] — [4.5]), <r 2 = cr([3.6] — [5.8]), and cr 3 = cr([4.5] — [8.0]) are combined 
errors, added in quadrature. 


I 11 the second phase, sources with 24 /irn data in the remaining pool are reexamined. 
Sources with colors of [5.8] — [24] > 2.5 or [4.5] — [24] > 2.5 are classified to be “transition 
disks”. With the potential contaminants and YSOs identified above excluded, there still 
remain some sources with bright 24 //in emission. For sources lacking valid photometric 
data in one or more of the IRAC bands, they are picked out as additional Class I type 
YSOs once fulfilling [24] > 7 and [X] — [24] > 4.5 mag, where [ X ] is the longest wavelength 
IRAC detection that we have. 


The classification scheme of 


Gutermuth et al. 


(120091 ) includes three phases, using 


2MASS data to identify more YSOs in addition to the aforementioned methods. However, 
in our case as G46 is much farther than the sources in 


Gutermuth et al. 


020 091 ). 2MASS 


photometry would be severely affected by foreground interstellar extinction and the use of 
these near-IR data to identify YSOs would induce heavy contamination from foreground 
field stars. Thus, we do not use the near-IR bands in our YSOs search. 


The above identification procedures have resulted in 106 YSOs. And they are classified 
into 22 Class I, 60 Class II, and 24 “transition disks” objects. A summary of the results 
is given in Table Qj We note that the color based classification scheme would lead to 
misidentihcation. A Class II YSO viewed at high inclination would show features resembling 
that of a Class I source. An edge-on Class I YSO can have similar infrared color of a 
Class 0 source. Thus, all Class I, II, and “transition disk” sources identified in this paper 
are YSO candidates. The distribution of these YSO candidates on distinct color-color 
spaces are presented in Figs. [ 8 ] and [9j It is important to note that we also use the point 
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sources extracted from 
identical results. 


Gutermuth &; Heverl (120151 ) to identify YSO candidates and reached 


3.3.2. Active Star Formation in GRSMC GO 46 .34-00.21 

Figure ITUT a) shows the large scale spatial distribution of the 106 YSOs identified in 
the region. Among the 82 Class I and Class II YSOs more than 40 (about 50 percent) 
are located in projection onto the molecular cloud GRSMC G046.34-00.21 (see Fig. UJ) 
which fills only the 20 percent of the whole survey held. Such a concentration of YSOs is 
suggestive of active star formation taking place in the cloud. By the other hand, from the 
24 “transition disk” found in the whole held, no one is detected in association with the 
molecular gas, indicating the relatively young nature of the YSOs embedded in this cloud. 
Interestingly all YSOs associated with the molecular cloud are placed in projection between 
the open border of G46 and the head of the pillars. Moreover, there are not young sources 
in the farthest part of the cloud behind the heads of these pillars (see Fig. flOTb)) . 

From Figure ITOV ad it can be appreciated only two concentrations of young objects over 
the whole held, one centered at 19 h 16 m 30fi, +11°51'00" (J2000), consisting of 15 Class II 
type sources and the other one centered at 19 h 17 m 00B, +11°47 , 00" (J2000) mostly composed 
by Class I type YSOs (9 sources), which are located just ahead the pillars-like features. A 
closeup view of the pillar-like features region is shown in Figure fTUlfb) . There is a point 
source revealed at 24 /irn which is too faint at IRAC bands and not listed in the GLIMPSE 
catalog. We have marked this 24 /irn point source using a red circle in Figure HOf bh Weak 
at short wavelengths, this point source would be younger than other ones and could be a 
Class 0 candidate. 

The fact that the Class II concentration is located closer to the open border of G46 
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than the Class I group strongly suggests an age gradient in the YSOs distribution. By 
the other hand, the absence of young sources inside and behind the pillar-like structures 
could show that the Hll region influence has not reached the molecular material behind 
the pillars, which is in agreement with the result of Section 13.2.11 We suggest a scenario 
where the propagation of ionization radiation escaping from the Hll region has triggered 
the star formation observed in the molecular cloud through RDI mechanism and then has 
stalled at the surface of the pillars heads. We can not discard that a growing density of the 


1BL reaches equilibrium wit 


i that of the cloud and thus, the shock front will continue their 


propagation into the cloud (Thompson et al 


2004a). 


Additionally, in the region surveyed with the JCMT (see yellow rectangle in Fig. [4]) 
there are 5 YSO candidates. One of them, a Class 1 source, lies exactly at the center of the 
dense molecular clump mapped with the HCO + and HCN J=4-3 (see Fig.lHJ), which is in 
agreement with the molecular gas conditions studied in Sec. l3.2.2l 


4. Conclusions 

We have analyzed the Hii region G46.5-0.2, its molecular environment and the young 
stellar objects placed in a wide region around it, searching for evidence of their physical 
connection and the possibility of induced star formation. We found a rich combination of 
mutual influences, underscoring the important role of Hll regions in favoring the formation 
of new stars not only on their immediate vicinity but also farther than a distance equivalent 
to its own radius. 

We characterized that the Hii region, located at about 4 kpc and with a radius of ~ 4.6 
pc, is probably excited by an 07V star. Its shape resembles a horse-shoe, with the bright, 
thick border towards the east-northeast and the open portion towards the southwest side. 
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When mapping the distribution of the molecular gas associated with G46 we found that the 
Hu region is located close to the edge of the GRSMC G046.34-00.21 molecular cloud and, 
curiously, G46 instead of opening in the direction of lower ambient density it does exactly 
in the opposite direction, towards the cloud. Filamentary structure in the molecular cloud 
is observed, particularly in the 12 CO J=3-2 and 13 CO J=3-2 maps. Towards the open 
border of G46 do not appear any considerable molecular emission, suggesting the presence 
of a pre-existing region with scarce molecular gas or that the UV-photons have carved the 
molecular cloud. Besides, close to the end of the observed filamentary structures in the 
cloud, the infrared images reveal the existence of pillar-like structures pointing towards 
the Hll region open border, which are associated with some narrow molecular filaments 
characterized using the 13 CO J=l-0 line. From a pressure balance study we found that 
the internal pressure of the neutral gas in the pillar-likc features heads are larger than 
the external pressure due to the ionized gas stalling at their tips, discarding that the RD1 
mechanism is on going. In addition, we analyzed a compact molecular clump located over 
the south border of G46, obtaining a mass of ~ 375 M 0 and a density of about 10 4 cm -3 , 
with a higher density towards its center, where the HCN J=4-3 line is detected. 

Taking into account that the expansion of the Hll region itself, and the injection of 
extra energy into the molecular cloud can drive turbulence and trigger star formation we 
searched for candidates of YSOs and classified them according to their evolutionary stage. 
We identified two main concentrations of young objects over the whole region, one, closer 
to the G46 open border consisting of Class If type sources, and other one mostly composed 
by Class 1 type YSOs located just ahead the pillars-like features strongly suggesting an age 
gradient in the YSOs distribution. 
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Fig. 1.— Two-color composite image of a large area towards the Hu region G46. In cyan 
the Spitzer -IRAC 8 /im emission, and in red the Spitzer -MIPSGAL emission at 24 /im . 
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Fig. 2.— Two-color composite image (8 /irn = cyan and 20 cm = red). The symbols represent 
the only eight point sources with JHK near-infrared and RB optical bands measurements 
found in the region (dashed circle). The diamonds (sources #6 and #8) indicate the location 
of the most likely candidates to be the exciting stars of G46. The beam of the 20 cm emission 
from MAGPIS is included in the top right corner. 
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Fig. 3.— Fitting of the spectral energy distribution for sources 7^6 and 7 ^ 8 . The black curves 
correspond to the best photospheric model of Kurucz (1979). The dots represent the data 
used for the fitting. 


12:00:00 


5 11:48:00 

a 



19:17:34.7 19:16:42.4 

Right ascension (J2000) 


Right ascension (J2000) 



Fig. 4.— Left: Spitzer -IRAC 8 pm emission with contours of the 13 CO J=l-0 emission 
integrated between 40 and 62 km s _1 , with levels of 10, 15, 20, and 25 K km s -1 . The yellow 
circle is the beam of the observations. The rms noise is about 1 K km s -1 . The rectangle 
represents the region studied with the JCMT data. Right: Integrated 13 CO J=l-0 from the 
GRS towards the pillar-like features observed in 8 pm. The white contours represent the 
13 CO emission integrated between 54.0 and 56.5 km s ^ 1 with levels of 4.5, 5.0, and 5.5 K 
km s _1 , while the yellow ones represent the integration between 56.5 and 59.0 km s -1 with 
levels of 4 and 5 K km s^ 1 . The beam of the molecular data is shown in the top-right corner. 
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Fig. 5.— 12 CO, 13 CO, and C 18 0 J=3-2 emission integrated between 40 and 60 km s~ 3 dis¬ 
played in contours over the Spitzer -IRAC 8 fin r emission. The contours levels are 45, 58, and 
72 K km s -1 ; 13, 20, and 27 K km s -1 ; and 3.3, 5.0, and 6.7 K km s _1 , for the 12 CO, 13 CO, 
and C 18 0, respectively. The rms noises are about 2.5, 0.8, and 0.5 K km s _1 , respectively. 
The beam of the observations is shown at the bottom left corner of the surveyed region. 



Fig. 6.— Spitzer-IKAC 8 pm emission with contours of the J=4-3 transition of the 
HCO + (left) and HCN (right) integrated between 40 and 65 km s _1 . The contours lev¬ 
els are 2.3, 3.3, 5.0, and 8.3 K km s _1 for the HCO + , and 1.2, 2.0, and 2.8 K km s _1 for 
the HCN. The rms noise is about 0.3 K km s _1 for both emissions. The yellow rectangle 
represents the same area as shown in previous figures. 
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Fig. 7.— Left: CO isotopes spectra. The rms noises of each spectrum are: 250, 400, and 
200 mK, respectively. Right: HCO + and HCN spectra. The rms noises are 33 and 116 mK, 
respectively. All spectra were obtained towards the center of the dense molecular clump 
related to the IRDC 046.424-0.237 and to the millimeter continuum source BGPS G046.427- 


00.237. 
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Fig. 8.— IR.AC color-color diagrams. Class I, Class II, and “transition disk” YSOs are 
marked using red triangles, blue plus symbols, and cyan diamonds, respectively. Mean 
errors of colors are pr esen ted in the bottom right corner of each panel. The definition of loci 


are adopted from 


Gutermuth et ah 


(120091 1. 
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Fig. 9.— IRAC-MIPS color-color diagrams. Class I, Class II, and “transition disk” YSOs 
are marked using red triangles, blue plus symbols, and cyan diamonds, respectively. Mean 
errors of colors are presented in the bottom right corner of each panel. 
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Fig. 10.— (a) Distribution of YSOs. Class I, Class II, and “Transition Disk” YSOs are 
marked using red triangles, blue plus symbols, and cyan diamonds, respectively. The 
grayscale shows emission at 8.0 /mi. (b) A closeup view of the pillar-like features region. 
The background presents emission at 24 /im. 
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Table 1: Exciting star(s) candidate(s) for G46. 


# 

2MASS Designation 

J 

H 

K 

B 

R 

A v [mag] 

T e//[K] 

2 

Xperpoint 

Spectral type 

3 

19172025+1156369 

12.508 

12.130 

11.942 

16.80 

14.70 

4-5 

~26000 

1.08 

later than B0 

6 

19172053+1154438 

12.498 

12.108 

11.9815 

15.70 

14.10 

3.5-4.5 

~30000 

1.28 

~09.5V 

8 

19172239+1155353 

10.094 

9.353 

8.914 

18.40 

14.80 

4-5 

~35000 

1.16 

~[07-07.5]V 


Table 2: Line parameters for the molecular lines towards the center of the southern molecular 
clump. 


Molecular line 

T mb peak 

(K) 

V LSR 

(km s -1 ) 

Av (FWHM) 

(km s' 1 ) 

12 CO (3-2) 

13.90i0.80 

52.50i0.05 

2.95i0.15 


1.93T0.33 

56.10i0.15 

0.95i0.25 


2.92T0.25 

58.70i0.30 

6.10i0.80 

13 CO (3-2) 

7.83T0.25 

52.60i0.01 

3.47i0.15 

C 18 0 (3-2) 

2.75T0.33 

51.80i0.15 

2.80i0.40 

HCO+ (4-3) 

2.16T0.15 

53.70i0.04 

3.15i0.10 

HCN (4-3) 

0.45i0.12 

53.10i0.40 

5.50il.00 


Table 3: Radex results from the HCO + and HCN J=4-3 lines using Tk = 20 K and the 
indicated nn 2 . 


input nn 2 

(cm -3 ) 

N(HCO+) 

(cm -2 ) 

^00+ 

N(HCN) 

(cm -2 ) 

ThCN 

A+CN/HCO+ 

10 5 

1.9 xlO 14 

7.50 

8.4 xlO 14 

4.50 

4.42 

10 6 

1.7 xlO 13 

0.70 

6.5 xlO 13 

0.37 

3.82 

10 7 

5.4 xlO 12 

0.14 

6.8 xlO 12 

0.09 

1.25 




Table 4: Source Counts in the YSO Search 


Sources 

G46 Field 

In GLIMPSE Catalog 

19210 

Valid in all Four IRAC Bands’ 1 

4491 

With 24 jiiii data 

338 

Rejected Contaminants 

883 

Class I Type 

22 

Class II Type 

60 

Transition Disks 

24 

Total YSO Candidates 

106 


Sources with photometric uncertainties no larger than 0.2 mag in all four IRAC bands. 



